Two-Qubit Gates for Resonant Exchange Qubits 
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A new approach to single-qubit operations using exchange interactions of single electrons in gate- 
defined quantum dots has recently been demonstrated: the resonant exchange qubit. We show that 
two-qubit operations, specifically the CPHASE gate, can be performed between resonant exchange 
qubits very straightforwardly, using a single exchange pulse. This is in marked contrast to the 
best known protocols for exchange qubits where such a gate requires many pulses so that leakage 
processes arising from the exchange interaction can be overcome. For resonant exchange qubits a 
simple two-qubit gate is possible because in this mode of operation energy conservation suppresses 
leakage. 



Since the seminal work of Loss and Divincenzo [T], 
the spins of individual electrons trapped in gate-defined 
quantum dots have been a very promising architecture for 
quantum computing [2HTT] . Various modes of operation 
have been implemented that pose different experimental 
challenges. An early theoretical advance was the reali- 
sation that if a single qubit is encoded in three electron 
spins then universal operations can be performed using 
modulated exchange coupling alone |12| . This in princi- 
ple removes the need for individually addressable electron 
spin resonance or large magnetic field gradients. There 
have been a series of experiments designed to pursue this 
avenue for quantum computing, despite the seeming in- 
efficiency of using three electron spins to encode a single 
qubit (HHII]- 

A key experimental challenge of the exchange qubit 
of [H] is the complicated pulse sequences that are used 
to implement two-qubit gates, and these have yet to be 
experimentally demonstrated. The pulse sequences are 
required to reverse spin flip transitions induced by ex- 
change coupling between qubits. DiVincenzo et al per- 
formed a numerical search to find an explicit pulse se- 
quence to implement a CNOT gate between the two 
qubits that required 19 exchange pulses in 13 time- 
steps [12]. This pulse sequence was subsequently shown 
to correspond to an exact analytic solution [TJ]. A vari- 
ant of this scheme, known as the decoherence free sub- 
system, has protection against uniform magnetic fields 
in all directions and recent results have found pulse se- 
quences for this encoding that are not significantly more 
complicated, requiring 22 pulses in 13 time-steps [l4] . 

Recently an alternative approach to exchange-only 
qubits has arisen, the resonant exchange qubit [151 I16j . 
As in the exchange qubit of [12] a single qubit is en- 
coded in the spin state of three electrons, each of which 
is trapped in an individual quantum dot. The qubit is 
operated with significant exchange interaction between 
the three dots and universal single-qubit operations are 
perfotmed by rf gate pulses that modulate this exchange 
interaction. It has been shown both theoretically [15] 
and experimentally [1_6J that this qubit has several advan- 



tages, including first-order insensitivity to charge fiuctua- 
tions and reduced leakage error due to nuclear field fiuc- 
tuations. This is a very promising demonstration of a 
qubit provided that one can couple them efficiently with 
low leakage. Taylor et al. show that this can be done 
through dipole interactions between the qubits [15] , while 
we suggest the alternative of using exchange coupling be- 
tween nearby qubits. 

In this paper we demonstrate that a further advantage 
of the resonant exchange qubit is that exchange-based 
two-qubit gates can be performed very simply. The key 
idea is that the leakage processes that arise from spin 
fiips do not conserve energy in the resonant exchange 
qubit. During an exchange pulse between two qubits 
these leakage processes are very strongly suppressed and 
two-qubit gates such as CPHASE can be performed in 
a single pulse in a single time-step. Moreover, leakage 
can be made arbitrarily small by reducing the strength 
of the exchange coupling relative to the resonant qubit's 
energy splitting. This mechanism contrasts with the con- 
ventional approach to exchange-only qubits where the 
leakage processes are reversed using long pulse sequences. 
The physics that leads to non-trivial two-qubit gates with 
low leakage using exchange coupling has also been pro- 
posed as the mechanism for two-qubit gates in the con- 
text of spin cluster qubits [TT] [T5] and some of our analy- 
sis is similar to the discussion in those references. More- 
over, the idea of energetically suppressing spin flips for 
quantum gates in quantum dot qubits specifically has 
previously arisen in the slightly different context of mag- 
netic field gradients and one- or two-spin qubits [T9H2T] . 
Exchange-based entangling gates between two resonant 
exchange qubits would appear to be feasible in the near 
future and in the rest of the paper we investigate this 
scheme in more detail. 

A single resonant exchange qubit is a triple-dot sys- 
tem operated deep in the (1, 1, 1) charge state with one 
electron in each dot. A large uniform magnetic field is 
applied in the z-direction such that the 8 electron spin 
states are Zeeman split according to the z-component of 
their total angular momentum m^. Tunnelling between 
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TABLE I: Energy eigenstates of the triple-dot system that 
constitutes a single resonant exchange qubit. S is the total 
angular momentum quantum number of the three electron 
spins and m^ is the z-component of the total angular mo- 
mentum, J^ is the exchange splitting for the qubit and B is 
the Zeeman splitting energy for a single spin. 
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FIG. 1: Diagram of three geometries for coupling resonant 
exchange qubits. Solid lines indicate the large exchange cou- 
pling Jz within a single qubit. Thin lines indicate the weaker 
exchange coupling Jc that is turned on between qubits to 
produce two-qubit gates, a) Linear geometry, b) Butterfly 
geometry, c) Rectangular geometry. 



adjacent dots is tuned by applied gate voltages resulting 
in two exchange splittings J;, coupling dot 1 to dot 2, 
and Jr, coupling dot 2 to dot 3. As explained in |15| and 
in |16j . whose notation we largely follow, the qubit oper- 
ating point has J/ = Jr- The qubit states |0) — 4=(| fti 
) + I ;n) - 2| m)) and |1) = ^(1 m> - l ;n» are 
eigenstates of the Hamiltonian describing the exchange 
interaction between electrons in each dot. They have 
rriz = 1/2 and S — l/'i,, where S is the total angular 
momentum quantum number. The qubit states are split 
in energy by J^ = (J/ -|- Jr)/2, with |0) being the lower 
energy state. In [TH] Jz/h is in the range 0.2-2 GHz. 

A third eigenstate |Q> - ^(| ttl) + I Itt) + I m)) 
has mz = 1/2 and therefore the same Zeeman energy 
as the two qubit states, but S = 3/2. The five remain- 
ing spin states have different Zeeman splittings and are 
given by |Q|) = | ttt) and |0_|), |g_), |0_), |1_) which 
are obtained from IQi), |Q), |0) and |1) respectively, by 
flipping all spins. The energies of these states of the sin- 
gle qubit system are indicated in Table IT] (In the rest of 
this paper we will take h = 1 so that the entries in the 
table can be regarded as either energies or angular fre- 
quencies.) Single qubit operations are performed using 
oscillatory exchange pulses [111 [TB] . 

We now consider a device involving two such qubits, 
shown schematically in Figure [Tl The first qubit involves 
dots 1,2 and 3, while the second qubit will involve dots 4,5 
and 6. We will label the first (second) qubit A (B) and 
use JzA (Jzb) to indicate the exchange splitting between 
the qubit states. We will not assume that JzA — JzB- 

Consider the system with no coupling between the 
qubits, the eigenstates are just the 8^ = 64 tensor prod- 
ucts of the states in Table|l| The four possible qubit states 
are in the m^ = 1 subspace of the six-spin system. This 
subspace has 15 states, so there are 11 leakage states with 
the same Zeeman energy that can in principle be coupled 



by inter-qubit exchange pulses. Table [II] indicates these 
eigenstates and energies. Whatever exchange coupling is 
turned on between the qubits, it will also conserve to- 
tal angular momentum. Six linear combinations of these 
leakage states have total angular momentum 5' = 2 or 
3 and so cannot couple to the qubit states which have 
S — \. Moreover exchange couplings will not mix states 
that have different values of ruz- As a result we may 
focus attention on the block of states with TOz — 1 and 
neglect other states from now on. 
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TABLE II: Energy level diagram for two resonant exchange 
qubits in the subspace with m^ = 1. J^j is the exchange 
splitting of qubit I — A,B. Without loss of generality we 
have assumed that Jzb > JzA- Note that no qubit state has 
the same energy as any leakage state. All energies are shifted 
by the Zeeman energy which is —B for all the states in the 
table. 

The most important feature of Table [ll| is that none of 
the leakage state energies correspond to qubit state ener- 
gies. So, for example, when JzA = JzB — Jz each qubit 
state is detuned from all leakage states by an energy of at 
least Jz/2. This means that spin flips that result from ex- 
change coupling between the qubits will be energetically 
unfavourable and leakage rates will be reduced so long 



as the exchange couphng between qubits remains smaller 
than the resonant qubit's energy splitting J^. This will 
mean that two-qubit gates will take a time long compared 
to 1/Jz- However, the single qubit operations for reso- 
nant exchange qubits must already run slowly compared 
to this timescale, so this is not a restriction in practice. 
We will show below how to perform non-trivial two-qubit 
gates in a time of order 25(27r/ Jz) for which leakage er- 
ror is not expected to be the dominant source of error. 
For the parameters of |16| this time is around 65 ns and 
compares favourably with the time taken for single-qubit 
operations. 

We will now explain several approaches to implement 
specific two-qubit gates using exchange coupling between 
the two sets of three dots. The qubits are coupled 
by inter-qubit Heisenberg exchange coupling Hij{Jc) = 
Jcio'xiCxj + o'yiO'yj + (Jzi<Jzj — -^)/4 of the ith dot with 
the jth dot, where i € {1, 2, 3} and j € {4, 5, 6} and J^ is 
the corresponding exchange energy and is controllable by 
adjusting some combination of gate voltages. (One could 
also turn on several such couplings.) 

We will consider three different geometries, as indi- 
cated in Figure [T] the linear, butterfly and rectangular 
geometries. In the linear arrangement one turns on an 
exchange interaction H34 between the third and fourth 
dot. In the butterfly arrangement one couples the second 
and fifth dots (-ff25)- Finally in the rectangular arrange- 
ment one can couple three sets of of neighbouring pairs 
of dots through H14 + H25 + H^q. These geometries have 
different advantages and disadvantages in terms of the 
simplicity with with which certain gates can be obtained 
and their susceptibility to noise. However, the choice of 
geometry will also depend on which arrangement of dots 
can most easily be tuned into the (1, 1, 1, 1, 1, 1) charge 
configuration and which geometry best allows for cou- 
pling large arrays of qubits. 

As noted above leakage processes will be suppressed 
when Jc <C Jz so we first analyse these couplings in low- 
est order perturbation theory. This requires us to cal- 
culate the level shifts of the two non-degenerate qubit 
states 1 00) and |11), and to find matrix elements of the 
coupling Hij on the subspace spanned by the two pos- 
sibly degenerate states |01),|10). To write the answer 
explicitly we restrict attention to the qubit subspace and 
define logical Pauli operators such that crzA\0)A = \Q}a 
and azA\i)A = -|1)a, crxA\0)A = |1)a, etc, and likewise 
for qubit B. Then the Hamiltonian Hq describing the 
two uncoupled resonant exchange qubits can be read off 
from Table II and is as follows 
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The overall level shift given by the first term is not im- 
portant so long as we can ignore the leakage states. For 
each of the dot geometries the coupling of the two qubits 
in lowest order perturbation theory can be summarised 



by the following effective Hamiltonian 

He = 5Jo+5Jz{(yzA+OzB)/'2 

+ JzzCrzA'^zB + J±i(^xA<7xA + CTyAC^yA)- (2) 

where the coefficients SjQ,dJz,Jzz,J± are each propor- 
tional to Jc, the applied exchange coupling, and are 
straightforward to calculate. Table [Hi] gives explicit val- 
ues for SJz , Jzz ; "^j^ in each of the three geometries at 
lowest order in Jc. The first term describes a level shift 
of the qubit states relative to the leakage states. The 
second term describes a shift of the qubit level splittings 
JzA,JzB- This will generate single qubit unitaries that 
can be corrected if necessary. The final two terms cou- 
ple the two qubits and can be used to implement two- 
qubit gates. Beyond lowest order in perturbation theory 
the structure of equation pi will be preserved but the 
parameters have a more complicated dependence on Jc- 
Note that, similar analyses of Heisenberg interactions of 
clusters of spins appear in many places in the literature, 
for example [13 HI [52] . 
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TABLE III: Qubit coupling parameters arising from lowest 
order perturbation theory in each of the three geometries of 
Figure[l| Calculations assume that JzA — JzB- When \JzB — 
Jza\ 3> Jc we find J± — for all geometries, but the other 
entries in the table are unaffected. 

We will now describe how to use the interaction 
of equation (pi) to implement specific gates. We will 
largely restrict our attention to CPHASE gates, although 
other choices are possible. CPHASE is the unitary 
diag(l, I, I, — I). In general we will directly obtain a 
gate that is equivalent to CPHASE up to local unitaries, 
(CNOT is one such gate), but since universal single qubit 
operations are known for the resonant exchange qubit 
this is enough to provide a universal gate set. The first 
general approach, the DC scheme, involves slowly turning 
on and off a small Jc- With an interaction time of order 
1/Jc we obtain a CPHASE gate. We will also consider 
AC schemes where Jc{t) is modulated at high frequency. 

The butterfly configuration provides an especially sim- 
ple CPHASE due to the enhanced symmetry in this case. 
For this configuration we find Jzz — Jc/^ and Jl = so 
the two qubits have an effective Ising interaction and an 
exchange pulse Jc{t) having area (j) = J Jc{t)dt = 97r/4 
results in a unitary U{(j)) that is equivalent under local 
unitaries to CPHASE. 

For the linear and rectangular geometries the most 
straightforward way to obtain the CPHASE gate requires 
that the two qubits have rather different exchange split- 
tings \JzB — Jza\ ^ Jc, in which case the qubit states |0f) 



and 1 10) are no longer degenerate and so Ji = 0. Again 
a CPHASE results from a single pulse Jc{t) having the 
correct area such that / Jzz{t)dt = tt/4. If it is preferable 
to operate in the regime where Jza = JzB then achiev- 
ing a CPHASE gate will require some single-qubit gates 
to echo out the effect of the J±_ term in equation (pi). 
Methods to find explicit pulse sequences exist, see for ex- 
ample |23j and references therein. One simple example 
is given in the Supplementary Material. 

There are three main sources of errors for this gate: 
charge noise (which leads to noise on J^), magnetic field 
noise, and leakage. Whereas single qubit operations for 
the resonant exchange qubit are first-order insensitive to 
charge noise, it is clear from Equation ^ that fluctu- 
ations in Jc will lead to correlated dephasing of the in- 
dividual qubits through the SJz oc Jc term as well as 
to fluctuations in the pulse area </>, reducing the fidelity 
of the gate. Specifically, low-frequency charge fluctua- 
tions will be significant as in the singlet-triplet qubit [8] . 
Note that since the rectangular geometry has SJz = it 
may perform better with regard to charge noise than the 
other geometries. Nevertheless, experimental experience 
in two- and three-dot qubit devices |71 |S1 IH] shows that 
high fidelity exchange pulses of the required length can 
be carried out in practice. Like single qubit operations 
of the resonant exchange qubit, the two-qubit gate will 
be affected by piezoelectric coupling to phonons [TS] and 
magnetic field fluctuations, specifically to low frequency 
fluctuations of the z-component of the Overhauser field 
gradient [Hj. These lead to independent damping and 
dephasing of the two qubits and are largely unaffected 
by the operation of the two-qubit gate. Finally, although 
the lowest order perturbation theory analysis given above 
does not have any leakage, this will occur in practice and 
the exact behaviour will depend on the pulse-shape cho- 
sen for Jc{t). Sufficiently smooth pulses with sufficiently 
small values of dJc/dt will greatly suppress leakage rela- 
tive to a square pulse. 

The error in the gate due to leakage is unlike the other 
sources of error in that it may be reduced by performing 
the gate over a longer time, using a smaller value of Jc- 
For example, we can choose to measure the leakage er- 
ror by maximising the probability of being in a leakage 
state after the gate over all input qubit states to find the 
maximum leakage probability pl ■ If there were no other 
source of error the worst case fidelity of the gate would 
then he Fl = 1 — PL (This is a lower bound on the gate 
fidelity.) Simple numerical calculations for the butter- 
fly geometry, which has the lowest leakage error of the 
three geometries, show that when Jc < O.ISJ^ the worst 
case leakage probability pl is less than 1%, even with a 
square pulse. With Jz — 2tt x 0.36 GHz as in [16] this 
leads to a gate time to — Qtt/AJc ~ 21 ns. This compares 
favourably with the time required for single qubit oper- 
ations in [TB]. Leakage errors reduce rapidly as the gate 
time is increased. So for square pulses with Jc < 0.05 Jz 



having a gate time of around 63 ns we already find from 
numerics that the leakage error is below 0.1% . This 
shows that with practical gate times leakage error may 
be reduced to the point where it is expected to be much 
less significant than other sources of error. The Supple- 
mentary Material has further discussion of leakage. 

We briefly consider an AC scheme, where additional 
control can be obtained by choosing an rf modulation 
of the the exchange coupling Jc- For simplicity, we will 
consider the rectangular geometry only and choose the 
two qubits to have different exchange splittings so that 
AJz — JzB — JzA is larger than the coupling Jc- We 
can then choose Jc{t) = Jcoit) + Jc\{i) cos{A Jzt) where 
Jco and JcA are both slowly varying functions of time. 
Note that we must have \Jca\ < JcO since the exchange 
coupling is always positive. Substituting this into Equa- 
tion (pl) and making the usual rotating-wave approxima- 
tion, which is valid when Jco, | Jca| ^ AJ^, we find the 
following interaction picture Hamiltonian 
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-crzAf^zA 



JcA 

~2A 



{f^xAf^xB + CTyACTyB)- (3) 



Since Jco and JcA can be independently controlled this 
interaction allows greater flexibility than in the DC 
scheme where the parameters in equation (pl) are all pro- 
portional to a single control parameter. So this interac- 
tion results in a CPHASE gate in a single pulse of area 
(j) — Qn/4: with Jca — and a SWAP gate in two pulses by 
choosing Jca ~ JcO and = 67r followed by a CPHASE 
gate with JcA = 0. 

Since we cannot set Jco — without also having 
JcA = 0, gates based on this scheme remain sensitive to 
low frequency charge noise (fluctuations in J^o) as well 
as high frequency noise (fluctuations of Jca)- We are not 
aware of an exchange-based two-qubit coupling for the 
resonant exchange qubit that is intrinsically insensitive 
to noise in Jc- However, in the AC scheme echo pulses 
are available that achieve insensitivity to such low fre- 
quency charge fluctuations, an example is given in the 
Supplementary Material. 

In summary we have described a number of protocols 
for performing two-qubit gates on the recently proposed 
resonant exchange qubit. For the most convenient pa- 
rameters and device geometries we have shown how to 
implement CPHASE in a single exchange pulse. We leave 
quantitative estimation of gate fidelities in the presence 
of realistic noise to future work [21] . 
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SUPPLEMENTARY INFORMATION 



Echo Pulse Sequences 



In the main text we claimed that the model for two-qubit interaction in the DC scheme given by 

He = SJo + 5J^{a2A+crzB)/2 



(4) 



allowed for echo pulses to produce CPHASE gates even when J_l 7^ 0. 

The simplest example is in the rectangular geometry where 5Jz = and all terms in the Hamiltonian commute. If J±^ 
were zero we could implement CPHASE with a pulse having (po = J Jc{t)dt = 67r/4. The J± term in equation (|4| may 
be echoed out by breaking the Jc pulse into equal two pieces, implementing the pulse sequence azAU{(j)o/2)azAUX4'o/'^)- 
The single qubit gates azA can be implemented using suitable Rabi pulses on the first qubit. Due to these single qubit 
echo pulses, the phase acquired as a result of the J± term in the first coupling pulse is unwound in the second pulse. 
This four stage sequence results in a CPHASE. In the linear arrangement of dots both 6Jz and J± are non-zero and 
more complicated pulse sequences are required to obtain a CPHASE gate. 

In the main text we furthermore claimed that the model for two-qubit interaction in the AC scheme and the 
rectangular geometry given by 



Hint — 



6 



cO 



-CTzACTzA 



JcA 

~2A 



{cTxACrxB + CTyACTyB), 



(5) 



allowed for echo pulses to produce CPHASE gates while echoing out the, possibly noisy, phase due to Jco- 

For example, we can choose a pulse with JcA — JcO and pulse area — J JcA{t)dt — Stt. The gate sequence 
CxAU{4>)axAU{(j)) results in a gate that is equivalent to CPHASE and is insensitive to low-frequency charge fluctua- 
tions. The single-qubit gates ct^a can be implemented as a single Rabi 7r-pulse on the first qubit. Due to these single 



qubit echo pulses, the phases acquired as a result of the Jco term and the cfyACyB term in ([5]) during the first coupling 
pulse are unwound in the second pulse. The coupling pulses here, U{(f)), are each about twice as long as directly 
implementing CPHASE in the DC scheme. Assuming conservatively that each single qubit pulse takes around the 
same time as such a direct CPHASE, this gate sequence is about 6 times longer than the direct implementation. 



Symmetry analysis of coupling 

We present a more detailed discussion of the symmetries of the coupling Hamiltonian in the butterfly geometry in 
order to justify the statements in the main text about leakage. The extra symmetry in this geometry both acts to 
reduce leakage errors and enables more detailed analytical calculations. 

As discussed in the text this geometry involves the coupling Hamiltonian H25{Jc)- Notice that this coupling is 
invariant under swapping the outer dots of either qubit, that is we could swap 1 and 3 or 4 and 6 and have no effect on 
the Hamiltonian since it does not act on any of those qubits. These symmetries are also symmetries of the uncoupled 
two-qubit system. The coupling is also symmetric under swapping the two qubits. This symmetry under swapping 
qubits is only a symmetry of the full model when J^a = JzB ■ 

We should organise the eigenstates of the two-qubit system, as far as possible, according to their parity under 
these swap operations. If we consider just a single qubit we can notice that the qubit Hamiltonian is invariant under 
swapping dots 1 and 3 and that |0) = ^(| tti> + I ltt> - 2| tit)) and |Q) = ^(| ttl> + I Itt) + I tit)) have even 
parity under this swap while |1) = -^{\ fti) ^ | Itt)) has odd parity. [SW APi:i\l) = -|1)). Parities of the different 
tensor products of these states are straightforwardly determined. 

As discussed in the text all of the exchange couplings are invariant under global spin rotations and thus we 
are guaranteed that energy eigenstates can be chosen to have definite total angular momentum quantum numbers. 
Although this is not true of all the tensor product states discussed in the text it is not difficult to find the correct 
linear combinations of states with the same energy. We define the following energy eigenstates which all have definite 
total angular momentum 



\E) 

\E3) 
\EA) 

\F) 
\F2) 

\G) 
\G2) 

\K) 
\K2) 

\L) 
\L2) 




These eigenstates have the foUowing symmetry properties 



State 


Energy + B 


Parity 


AngularMomentum S 


Swap 


\E) 





++ 


1 


+ 


\E3) 





++ 


3 


+ 


\EA) 





++ 


2 


- 


\F) 


-Jza/2 


-+ 


1 




\F2) 


-Jza/2 


-+ 


2 




\G) 


-Jzb/2 


+- 


1 




\G2) 


-Jzb/2 


+- 


2 




|1,1) 


-{JzA + Jzb)/2 





1 


+ 


\M) 


-^Jza/2 


++ 


1 




\M2) 


-'■iJzAl^ 


++ 


2 




\N) 


-iJzB/2 


++ 


1 




\N2) 


-iJzBl2 


++ 


2 




|o,i) 


-{iJzA + JzB)/2 


+- 


1 




|i,o) 


-{JzA + iJzB)/2 


-+ 


1 




|0,0) 


-HJzA + Jzb)/2 


++ 


1 


+ 



The parity column indicates the parity under swapping the first and third dot and also the fourth and sixth dot as 
discussed above. Swap indicates the parity of the state under swapping qubits where this is well defined. 

Leakage 

In this section we justify the statements about leakage made in the main paper. We consider just the butterfly 
geometry. 

The exchange gate coupling Hamiltonian H25 can only couple eigenstates with the same total angular momentum 
and parity quantum numbers. So for example {E3\H25\ij) = where \ij) is any qubit state, because S = 1 for each 
qubit state and \ES) has 5 = 3. On the basis of their total angular momentum we conclude that the six states 
\E3), \EA), \F2), \G2), \M2), \N2) do not couple to the qubit states during an exchange pulse. The existence of these 
uncoupled leakage states was mentioned in the main paper. The five relevant leakage states are \E), \F), \G)\M)\N) 

Each of the qubit states has a different parity. As a result we find for qubit states {ij\H25\kl) = whenever 
ihj) 7^ (^1 ^nd we can conclude that when restricted to the qubit subspace H25 has no off-diagonal matrix elements. 

The state |11) is the only state with parity and so it doesn't couple to leakage states at all. Since \G) is the only 

leakage state with parity H — and 5 = 1 it is the only state coupled to |0, 1) by the exchange pulse. When J^a = JzB 
the energy difference between these states is iJz/2. Likewise |1,0) couples to F and the energy difference between 
these states is also 2>Jzl2. Finally |0,0) couples to the three remaining leakage states |_E), |M), |7V) that are different 
in energy by SJ^, iJz/2 and 37^/2 respectively. 

So in the butterfly geometry the enhanced symmetry means that each computational basis state couples indepen- 
dently to an orthogonal leakage subspace. As a result we may write a completely general solution for the mapping of 
the computational basis states under any exchange pulse Jc{t) as follows 



|0,0) - 


-> Vl - Pooe-^'^'"' |00) + V^ILoo 


|o,i> - 


-> ^Jl-poie-"^"'m + VP^i\Loi 


|i,o> - 


-> ^Jl-pwe-"^'°m + ^pra\Lio 


|0,0) - 


-> e-^'^^HOO) 



Where \Lqq), |Loi) and \Liq) are orthogonal leakage states. The leakage probability associated with the computational 
basis state \ij) is pij and because of the orthogonality of the \Lij) the largest of these is the maximum leakage 
probability over all possible input two-qubit states. This is the maximum leakage error p^ discussed in the main text. 
The symmetry of the problem requires that poi = Pio £^nd typically we find that poi > Pm, so usually p^ = Poi- 




FIG. 2: Plot of minimum gate fidelity Fl and worst case leakage probability Pl for a two-qubit gate having cj> = J Jcit) — 97r/4 
and JzA = JzB in the butterfly geometry, using a square pulse of time to- Time is normalised by tq = 2-71 /J^, the Bohr period of 
each resonant exchange qubit, which is around 3 ns for the parameters of [16j . The dashed horizontal line indicates a minimum 
fidelity of 0.999 or 0.1% error. 



To obtain a gate equivalent to CPHASE up to local unitaries one chooses the pulse area so that exp[— z(</joo + <(5ii — 
Vol ~ ^w)] = ^1- Once again due to the orthogonality of the leakage states, the minimum fidelity of this gate over all 
possible input two-qubit states is F^ = 1 — p^ in the absence of other sources of error. These quantities are plotted 
as a result of numerical integration for a square pulse of fixed area in Figure 1. These results are the basis of the 
quantitative statements about leakage in the main text. As noted in the main text fidelity is greatly improved, for 
the same gate time, by using smooth pulses rather than square pulses |21j . 



